25-Kicked Quantum Rotors: Localization and 'Critical' Statistics 
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The quantum dynamics of atoms subjected to pairs of closely-spaced 5-kicks from optical poten- 
tials are shown to be quite different from the well-known paradigm of quantum chaos, the singly-5- 
kicked system. We find the unitary matrix has a new oscillating band structure corresponding to 
a cellular structure of phase-space and observe a spectral signature of a localization-delocalization 
transition from one cell to several. We find that the eigenstates have localization lengths which 
scale with a fractional power L ~ ft - ' 75 and obtain a regime of near-linear spectral variances which 
approximate the 'critical statistics' relation Y,2{L) ~ xL ~ |(1 — v)L, where v ss 0.75 is related to 
the fractal classical phase-space structure. The origin of the v w 0.75 exponent is analyzed. 
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The 5-kicked quantum rotor (QKR) is one of the most 
studied paradigms of quantum chaos. Its implementation 
using cold atoms in optical lattices provided a con- 
vincing demonstration of a range of quantum chaos phe- 
nomena including Dynamical Localization 0], the quan- 
tum suppression of chaotic diffusion. Recently, an ex- 
perimental study [3| of cesium atoms subjected to pairs 
of 5-kicks (25-KR) showed surprisingly different behav- 
ior. The classical phase-space is chaotic but is made up 
fast diffusing regions which are partly separated by slow- 
diffusing 'trapping regions', where the classical trajec- 
tories stick; the classical analysis revealed a regime of 
anomalous diffusion corresponding to long-lived correla- 
tions between kicks. 

Further details are given in Q but we show here that 
the 25-KR has some unexpected quantum properties. We 
show that there is a cellular phase-space structure which 
arises from a novel oscillatory band structure of the cor- 
responding unitary matrix. One consequence is a new 
type of localization-delocalization transition not seen in 
the QKR, where states delocalize from single to multiple- 
cell occupancy; we show it has a clear spectral signa- 
ture. We have also found scaling behavior of the local- 
ization lengths associated with a fractional exponent, i.e. 
L ~ fir - 75 , whereas for the well-studied QKR, L ~ fir 1 . 
A similar exponent is found for the decay of return prob- 
abilities in the trapping regions, P(t) ~ t~°" 5 . We argue 
that the exponent 0.75 corresponds closely to the value 
obtained for the dominant exponent of the golden ratio 
cantorus |5j, |l| . We show that the spectral fluctuations 
(both the nearest-neighbor statistics (NNS) and spectral 
variances) show important differences with the QKR in 
regimes where the derealization of eigenstates is hin- 
dered by cantori bordering the cells. We find a regime 
approximating the form found in 'critical statistics': the 
number variances of the spectra are linear T,2(L) ~ %L 
for L » 1, where \ — 1/2(1 — v) < 1 and v ~ 0.75. 

The term 'critical statistics' arose originally in rela- 
tion to the Metal Insulator Transition (MIT) in systems 



with disorder |7J, |8| . A new universal form of the distri- 
bution of nearest-neighbor eigenvalue spacings, termed 
'semi-Poisson', Pis) ~ sexp— 2s was associated with 
the MIT 0- For critical statistics a very interesting 
connection has been established between the multifrac- 
tal characteristics of the wavefunctions and those of the 
spectral fluctuations H: the number variances of the 
spectra are linear Yj-JJS) ~ %L f° r i » 1- The slope, 
X — 1/2(1 — D2/D) < 1, was shown to be related to 
a fractal dimension D2 obtained from the second mo- 
ment of the wavefunction and to D, the spatial dimen- 
sion of the system. For integrable dynamics, in contrast, 
E 2 (L) = L while for a GOE, E 2 (L) ~ Ln(L). There 
is much current interest in so-called 'critical' statistics 
in non-KAM billiards (typically systems where the dy- 
namics would be integrable were it not for a disconti- 
nuity in the potential) [9|, which show multifractal scal- 
ings and linear variances related to Z?2- Below we apply 
the term 'critical statistics' in this broader sense, rather 
than the MIT critical point. Multifractal behavior has 
been demonstrated for Cantor spectra ^3] where the level 
density itself is not smooth. However, until now, criti- 
cal statistics have not been seen - and were not thought 
to be relevant to - KAM systems. These are systems, 
ubiquitous in many areas of physics, where the transi- 
tion to chaos as a perturbing parameter is increased is 
quite gradual. 

The Hamiltonian of the 26-KR is H(x,p) = ^ + 
K cos x J2 n $ — + — n T + e ); there is a short 
time interval e between the kicks in each pair and a much 
longer time ~ T between the pairs themselves. In exper- 
iments e ~ 0.01 - 0.1 < T and fi ~ 2 — > 1/4 in the usual 
rc-scalcd units 3] . 

A study of the spectral fluctuations of a time-periodic 
system involves a study of the eigenstates and eigenval- 
ues of the one-period time-evolution operator U (T, 0) . 
For the QKR, the matrix representation, in an angular 
momentum basis | I), has elements [7j m = Llf 7ee . Uj^ k = 
exp — il 2 Th/2 . Ji- m (jh). The 'kick' terms, J;_ m (-^) are 
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Bessel functions and give the matrix the banded form il- 
lustrated in Fig^a). Since Ji_ m (x) ~ if \l — m\ > x 
we define the usual bandwidth b = K/H. The resulting 
statistics are approximated by those of Band Random 
Matrix Theory (BRMT) [n| rather than of RMT: i.e. if 
the dimension of the U(T, 0) matrix is N to t, the statis- 
tics are Poissonian for Ntot 3> b; the eigenstates of the 
BRM are exponentially localized in I, with a localization 
length in momentum (p — IK) which equals L p ~ K 2 /h, 
so states separated in p by 3> L p will be largely uncorre- 
cted. 

For the 2£-KR, the corresponding matrix elements are: 



U, l: 



= e- il ^ T ~ 



(1) 

As the U(T, 0) matrix is quite insensitive to (T — e) the 
quantum dynamics largely depends only on two scaled 
parameters K e — Ke and h e = he, rather than on K, e 
and h independently: the remainder of the matrix is in- 
variant if K e and h e are kept constant. 

An analytical form for the bandwidth of U(T, 0) was 
obtained in : the bandwidth oscillates sinusoidally be- 
tween a maximum value b max — 2K/H for angular mo- 
menta I ~ 2nn/(he) and a minimum value b m in ~ for 
I ~ (2n+l)ir / (he) . These minima correspond to the trap- 
ping momenta p = lh ~ (2n + l)7r/e seen in experiment. 
The corresponding band-structure of U is illustrated in 
Fig^b): the band oscillates, and U is approximately 
partitioned into sub-matrices of dimension TV = '=% cor- 
responding to separate momentum cells. 

The key to our work is our ability to vary the transport 
between the cells (by opening/closing the classical fractal 
'gates' between them) separately from the degree of filling 
of each individual cell. We begin by introducing a 'filling 
factor' R where 



R = 



K 2 



Nh 2 2-kK 



(2) 



measures the degree of filling of a cell by a typical state in 
the absence of confinement. Clearly, if there is no trans- 
port between cells the states simply fill the cell uniformly 




FIG. 1: Left: Structure of time-evolution matrix U(T, 0), for 
the Quantum Kicked Rotor (QKR), in a basis of momen- 
tum states, exemplifying the typical band-structure of a Band 
Random Matrix (BRM). Right: U(T,0) for our system, the 
25-KR, showing the new form with oscillating bandwidth. Be- 
fore derealization, eigenstates are confined within a single 
'momentum cell' of dimension N. 



and even if R 3> 1, the localization lengths L p ~ Nh. 
We begin by defining a localized limit, where R <C 1 and 
L p <C Nh ~ — ; typical states are insensitive to the 
boundary conditions of the cell and this limit is Pois- 
sonian. At the other extreme, if we allow strong cou- 
pling between cells, we move to an opposite limit as an 
increasing proportion of eigenstates become delocalized 
over several cells. 

We now investigate the transport. A classical analysis 
0, 0] shows that if we take Ke <C 1 and expand initial 
momenta of the j — th trajectory about the trapping val- 
ues Pj = (2n+ l)n/e + Spj, we can show that much of the 
trapping region is given by a classical map quite similar 
to the well-known Standard Map: 

Pj+2 — Pj ~ N 2 — sin 2xj ~ K8p J e cos Xj (3) 

xj+2 ~ Xj +p j+2 T. (4) 

Over much of the trapping region, the second term 
in EqO is dominant, for parameter regimes of interest. 
Then the kick impulse has a ir/2 phase relative to the 
full map (the 2-kick map gives a pair of V'(x) — — sin 2 
type impulses) and a momentum dependent effective kick 
strength K' = KeSpj. A detailed study of classical 
phase-space shows that at low K' the resonance 
structure is locally quite similar to the Standard Map. 
Hence, in the regime KeSpj ~ 1 we expect the 'golden 
ratio' cantori, which result from the last invariant man- 
ifold of a KAM system, to provide the strongest barrier 
to transport 01 . Though we want Ke to be small, if 
Ke < 0.1 phase space becomes too regular. Hence here 
we find that the regime of interest is within the interval 
0.1 <K e < 0.7. 

We investigated the corresponding quantum transport 
by evolving a set of wavepackets <&(p, t) in time, (where 
$(p, t = 0) = S(p)) for a range of K e and h e , until the 
momentum spreading is arrested by dynamical localiza- 
tion at t ~ iff. The resulting probability distributions 
|$(p, t 3> tn)] 2 = N(p) have a characteristic 'staircase' 
structure, shown in Fig[2Ia). At each step there is a steep 
drop in probability: 

N(p) + = e- 2d N(p)_ (5) 

(where N(p)± represent probabilities before(-) and af- 
ter (+) the step) concentrated over the trapping region 
(~ 1/6 of a cell in every case Q). The staircase tracks 
an exponential envelope N(p) ~ exp — 2\p\/L exp , where 
L e xp — ~a • We average over several steps, to obtain d 
as a function of K e and h e . In Fig[2Jb) we show that, 
quite accurately, d oc hP c 75 j f(K e ) where f(K e ) is some 
function of the scaled kick-strength, K t . We estimate 

(6) 

This defines our transport parameter, d, and comple- 
ments the 'filling factor' R. The inner steps of the stair- 
case have been seen in the momentum distributions of 
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FIG. 2: (a) Final (t — > oo) momentum distributions, N(p), 
(slightly smoothed) for quantum wavepackets of the 25-QKR 
for K = 20, e = 0.0175 and ft = 1/8 and 1/30 respectively. 
N(p) for both the eigenstates and wavepackets shows a long- 
range 'staircase' form which on average follows the exponen- 
tial N(p) ~ exp — 2(p — p)/L exp where L exp /2 = 2n/he; the 
ft-dependence of L exp is determined by the drop in probabil- 
ity, d, at each step, (b) shows that Ln(2d) plotted against 
Ln(he) lies on straight lines of invariant K £ — Ke, with con- 
stant slope 0.75. Hence d oc (fie) ' 75 and L exp oc h~ ' 75 - 
in contrast to the well-known QKR result L exp oc ft -1 . The 
R = 1 border is shown: for R > 1 (below the line) eigenstates 
fill much of a single cell. The derealization border is shown 
at d ~ 2: for d > 2 (above the line), over 98% of the prob- 
ability of typical eigenstates is confined to a single cell; for 
d < 2, eigenstates begin to occupy multiple cells. Statistics 
are presented later in Fig. 3 for points corresponding to the 
dotted line. 

atoms in optical lattices so j though existing data 

is not in the critical regime, in principle the form of d is 
experimentally verifiable. 

We are unaware of another KAM system where a single 
power-law exponent is so dominant. Typically, power-law 
behavior is associated with mixed phase space behavior, 
where many competing exponents are found [To]. We 
note also that the value 0.75 coincides closely with one 
of the scaling exponents found in [j| for the golden ratio 
cantori: these were a « 0.65 (in the most unstable part 
of the cantori) and a w 0.76 (most stable regions). There 
have been previous studies of transport in a region near 



golden ratio cantori. These have found L ~ fr+°- e6 [l5j 
but only in a momentum band region 'local' to the can- 
tori. We note that the L ~ h +<J dependence is associated 
with the physical process termed 'retunnelling' [(| asso- 
ciated with cantori which are classically 'open' but for 
which h is too large to permit free quantum transport. 
An abrupt change to an L ~ h~ (T is observed when the 
cantori 'open' for quantum transport and L is determined 
instead by localization. It has been argued that the rea- 
son all previous studies have found L ~ fi+ - 66 6] is that 
the retunnelling transport favors the unstable direction. 
To our knowledge ours is the only example corresponding 
to the dominant exponent of the golden ratio cantori; we 
attribute this to the fact that we are always in a dynam- 
ical localization regime, and localization will select the 
most stable parts of the fractal cantori regions, where at 
low K, elliptic fixed points are found. 

We can now investigate the statistics as a function of 
the filling factor R and the inter-cell transport param- 
eter d. Full details are given in [4], but in brief: we 
considered two types of boundary conditions (BCs). (1) 
Periodic BCs, i.e. solving the problem on a 'torus' in mo- 
mentum space, a well-known procedure for the QKR 
(2) Open BCs, where we diagonalize U(T, 0) with N to t — 
10, 000, but N m 1000; we then assigned the i — th eigen- 
state to then — th cell if (2n+l)7r/e < (pi) < (2n + 3)7r/e 
and calculated statistics for each cell. In both cases we 
averaged over « 20 cells to improve significance. For pe- 
riodic BCs, eigenstates cannot escape from a single cell 
of width NH. For open BCs, however, they can delocalize 
onto neighboring cells. 

In Fig|2Ia), the ^(L) statistics are presented. These 
represent the variances in the spectral number density, 
T,2(L) = (L 2 ) — (L) 2 , where we consider a stretch of the 
spectrum with an average (L) levels. A fit to the best 
straight line in the range L ~ 5 — 40 yields an estimate 
of the slope x- I n FigEJb) we show the nearest neigh- 
bor P(S) statistics. We quantify the deviation of P(S) 
from Pp(S) and Pgoe(S), its Poisson and GOE limits 
respectively, with a parameter Q [l6| : 

1 Q- fo S °(P(S')-PGOE(S')dS>) 
J Sn (Pp(S')-P G oE(S'))dS'' 

Hence Q = indicates a Poisson distribution, while Q = 
1 signals a GOE distribution. We take So = 0.3. 

The results of FigOHdemonstrate that the statistics are 
not too sensitive to boundary conditions for d > 2 (see 
also 0|). However, for d < 2, while the states with peri- 
odic BCs (effectively restricted to a matrix of dimension 
N) move gradually to the GOE limit, for the open BCs, 
at derealization the statistics tend back to the Pois- 
son limit. This initially surprising behavior occurs be- 
cause, after derealization, one finds increasing numbers 
of states for which (pi) assigns them to the n — th cell, 
but for which much of the state's probability is actually 
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ties P(t) = \(ip(t = 0) |> 2 ~ t~ - 75 for wavepackets 
started in the trapping regions (see Q). The value of 
X ~ 0.125 corresponds to the value which would be ob- 
tained from the MIT relation, \ ~ 1/2(1 - D 2 ) with, 
if £>2 ~ 0.75; in the MIT, return probabilities with 
P(t) ~ t~° 2 were similarly found. We suggest that this 
represents a KAM analogue of behavior associated with 
'critical statistics'. 

In summary, we have shown that the behavior of the 
2<5-kicked system is rather different from the standard 
QKR. We have identified a spectral signature of the novel 
localization-delocalization transition that the system ex- 
hibits. We have also identified signatures of the frac- 
tal phase-space structure of the cell borders. The trap- 
ping regions may have applications in atom optics ex- 
periments, as a means of manipulating the momentum 
distribution of the atomic cloud. 

This work was supported by the EPSRC. We thank 
Shmuel Fishman and Antonio Garcia-Garcia for helpful 
comments and advice. 



FIG. 3: (a) shows the E 2 (L) statistics for the 25 - KR. 
Approximately linear variances, with slope \ ~ 0.125 are 
found in the intermediate regime. Solid line indicates the line 
^2(L) — 0.125L + c corresponding to the 'critical statistics' 
slope x ~ 1/2(1 — v) where v ~ 0.75 is the fractional expo- 
nent obtained previously, (b) shows corresponding nearest 
neighbor statistics for R = 0.8 and R = 1.3 of an intermedi- 
ate form, which is compared with the Semi-Poisson and GOE 
forms. The insets show the effect of the boundary conditions. 
Before derealization (d > 2) the results are insensitive to 
boundary conditions. At derealization, if the eigenstates are 
confined to a single cell (periodic BCs, points indicated by 
asterisks) the statistics make a transition to GOE; if they can 
delocalize (open BCs, points indicated by circles), the statis- 
tics tend to Poissonian behavior. The arrows indicate the 
d = 2 border. 



found in neighboring cells Hence states within each 
spectrum of N eigenvalues become progressively uncor- 
rected. This apparent failure of the procedure of assign- 
ing states to a given cell, in fact provides a rather good 
'marker' for the onset of derealization, and yields a clear 
'turning point' in both the NNS and ^(L) behavior. A 
detailed study of the multi-cell regime remains to be un- 
dertaken; models of the statistics for chaotic systems with 
non-uniform rate of exploration of phase-space |l7| may 
be relevant here. 

Of further interest here is an intermediate regime, 
found for both boundary conditions for R ss 1 and d > 2. 
Over a wide range of parameters and different cell sizes 
, we find approximately linear variances, for L> 1 and 
L < N, with slope x ~ 0.125. The inset of Fig^a) plots 
the values of x calculated along the vertical dotted line 
of Fij^Ib) for periodic and open boundary conditions. 
We note that a study of the decay of return probabili- 
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